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Abstract: We extend the approach of Banks, Myerson, and Kogut for the calcu- 
lation of the Wilson loop in lattice U{1) to the non-abelian SU{2) group. The original 
degrees of freedom of the theory are integrated out, new degrees of freedom are intro- 
duced in several steps. The centre group Z2 enters automatically through the appearance 
of a field strength tensor which takes on the values or 1 only. It obeys a linear 
field equation with the loop current as source. This equation implies that ff^^ is non 
vanishing on a two-dimensional surface bounded by the loop, and possibly on closed 
surfaces. The two-dimensional surfaces have a natural interpretation as strings moving 
in euclidean time. In four dimensions we recover the dual Abrikosov string of a type 
II superconductor, i.e. an electric string encircled by a magnetic current. In contrast 
to other types of monopoles found in the literature, the monopoles and the associated 
magnetic currents are present in every configuration. With some plausible, though not 
generally conclusive, arguments we are directly led to the area law for large loops. 
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1 Introduction 



It is now widely accepted, that confinement is due to the formation of a color electric 
string, and that magnetic monopoles play an essential role in this context. Up to now 
there is a lively activity in this field, illuminating the phenomenon from various sides. 
A particularly appealing approach is the one by Banks, Myerson, and Kogut [|l| which 
is now more than 20 years old. They considered the partition function for an electric 
current loop and derived step by step the appearance of monopoles by integrating out 
the original degrees of freedom and introducing new ones. The possibility to do this was 
restricted to the lattice U(l) model with the Villain action and some other simple 
models. The authors also remark, that the techniques do not generalize simply to non- 
abelian theories. 

We will start by applying and generalizing the methods of ref [ij. In a first step 
the SU{2) matrices on the links are explicitly parametrized by three angles ip,'d,ip. An 
appropriate decomposition of the SU{2) matrices allows the calculation of the trace in 
the plaquette action. After an expansion of exponentials into modified Bessel functions 
the integrations over the link angles can be performed. They lead to constraints for 
the new variables which were introduced in the expansions. Most of these variables are 
irrelevant and the summations can be performed after a suitable transformation. 

After this has been done, we are left with several integer variables which are restricted 
by constraint equations. The most important one is a Z2 field strength tensor /^j,. It 
lives on plaquettes and is either or 1. This tensor obeys an inhomogeneous linear 
field equation with the loop current as source. The solutions of the equation have a 
simple geometrical interpretation. The tensor /^^^ is non vanishing on a two-dimensional 
surface bounded by the loop, and possibly on closed two-dimensional surfaces. These 
surfaces have a natural interpretation as strings moving in euclidean time. There is a 
string which connects the two charges associated with the loop, and possibly a number 
of additional closed strings. The situation is particularly transparent for planar loops, 
where the layer on the minimal surface, corresponding to the straight string, plays a 
special role. For large loops we can use some general arguments and reasonings from 
statistical mechanics, like additivity of the free energy, and obtain the area law. The 
subtle question, whether a finite string tension survives in the thermodynamic limit and 
in the continuum limit, needs additional investigations. 

Our approach is purely analytical and non-perturbative, no gauge fixing was per- 
formed, and (3 was kept arbitrary. Lengthy calculations were done with the help of 
Mathematica. Nowhere any physical picture of what we expect was put in. It is the 
formalism itself which automatically leads to the appearance of a plaquette variable /^i, 
which is naturally associated with the world sheet of a string. If one could perform the 
integrations and summations over the remaining parameters one would obtain the ex- 
plicit string action. Even without doing this, the formalism clearly shows the appearance 
of the string picture and the origin of the confinement mechanism. 
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2 The partition function 



We are interested in the expectation value of a Wilson loop W, characterized by a closed 
current loop J: 

Z[J] = J Tr[WiJ)]exp[^Y.TrUMMU]. (2.1) 

Ppv 

The sum runs over all plaquettes Pp^, (with p < z/), while Upy{p) is the product of the 
four S'f/(2)-matrices on the links of the plaquette. We will perform rather extensive 
manipulations in the following, therefore we fix our notation here: 



p, q, r denote lattice points, 

fj.,^, p, X = 1, - ■ ■ ,d denote space directions, 

p ± yU is the lattice point next to p in positive or negative /^-direction, 

Pp denotes the link connecting p with p + p, (2.2) 

Ppu with /i < z/ is the plaquette determined by the links p^ and pu, 

0^(p) = {ip , -t} , !f) p{p) denotes three angles which parametrize the SU(2) matrices 

on the link p^, indices a, b generally run over ip, i), ip, 

Up{p) = UiQpip)) is the SU{2) matrix on the link p^, 

Upuip) = Up{j))U^{jp + p)Up{j) + u)U'^{j)) is the plaquette variable. 

As a first step we have to choose a parametrization for the link variables Up{p) in order 
to be able to do the group integrations. We proceed similarly as in previous work 
which applied the optimized ^-expansion on the lattice 0, 0- In our case the Euler 
parametrization is the most appropriate, i.e. we use 

U = e'^'^^e^^'^e'''^'^. (2.3) 

The Haar measure is proportional to sin(2^9). A possible choice of the parameter space 
is -vr <■?/'< vr, 0<v9<7r, Q < 'd < 7r/2. 

For the following it is convenient to extend this region. All integrals which appear 
contain functions of the Upy{p) which are periodic in the Euler angels. So we may use 
some symmetry relations which are easily seen from the decomposition of U into 1 and 
the am- The shift ip^ip — irjip^ip — n leaves U invariant. Therefore one can extend 
the 93-integration into the interval from — tt to tt, thus covering the group manifold twice. 
The further symmetry i) — > —d^ip ^ ip — 7^/2, ip ^ ip + tt/2 allows to extend the 'd- 
integration to the interval —tt/2 < 7? < 7r/2 if we continue the Haar measure as even 
function. Finally the symmetry d^d — ir^ip^ip — it allows the extension of the 
■{^-integration to the full interval. Therefore we take 

iJ(79) = ||sin(279)| (2.4) 

as Haar measure in the following, and use the common boundaries — vr < ip,{},ip < n. 
In fig. 1 we show the notation for the link variables of the plaquette Ppi^. 
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Fig.l. The plaquette Ppi, and the hnk variables. 
The parametrization according to (2.3) becomes 



Uu{p + p) 

utip) 



^iasipp{p) gj(T2i9p(p) gj(T3</Pp(p) 
^i(T3^v{p+p) g«o-2i9,^(p+p) ^i<T3'Pv{v+p) 
^-i(T3ipp{p+v) ^-icr2^p{p+u) ^-ia3ipp{p+u) 

g-j(T3(^^(p) g-j(T2i?y(p) ^-ia3ip^{p) 



(2.5) 



An appropriate technique for the further procedure, which was also extensively used in 
0] and 1^ , is the splitting of the matrix exponentials into sums of ordinary exponentials 
times projection operators, in general 



E 



PAm], with PJm) 



1 



l + sam)- 



(2.6) 



=±1 



The link variables U are parametrized as in (2.5), the plaquette variable Upu{p) therefore 
contains 12 factors. To every factor we apply the decomposition (2.6). This means that 
we need 12 summation indices s[B] = ±1, associated with the twelve angles in (2.5). At 
the end TrUp^ip) becomes a product of 12 factors, each being a sum of two terms of the 
type (2.6). So in total we have a sum of 2^^ terms. Each term of the sum consists of a 
trace T of a product of 12 projectors Ps{m), multiplied by an exponential. 

At the four corners of the plaquette one has a product of two projectors Ps{3)Ps'{3). 
They both project on the same subspace, therefore we get a non-vanishing result only 
if the neighboring parameters s and s' agree. So we must have s[ipp{p)] = s[ipi,{p + p)], 
s[ipu{p + p)] = s[ipp{p + u)],s[ipp{p + z/)] = s[ipi,{p)],s[ipu{p)] = s[ijp{p)], there are in fact 
not 12, but only 8 independent parameters. 

We enumerate the remaining 8 independent parameters s[G] in a consecutive way, 
starting with s['dp{p)]. In (2.7) we give the parameters Si, together with the angles to 
which they refer. According to the remarks above, the Si with even i belong to two 
angles. 
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S3 ^ -duip + p), S4<^ (Pu{p + p),Vp{p + Jy), 

S5 ^ '&p{P + ^): Sq^ 'll)p{p + v),ip^{p), (2.7) 
S7 ^ 1^u{p), S8<^ 1pu{p),1pp{p)- 

The number of terms in the expansion of TrUpiy{p) has now been reduced from 2^^ to 2^. 
We enumerate these terms by an index n, the order in which this is done needs not to be 
specified yet. Thus each Sj becomes a function of n, we denote it by Sin- The traces also 
depend upon n, we denote them by r„. Finally one can write TrUpi,{p) in the following 
form: 

TTUp,{p)^TT[Up{p)U,{p + p)U;{p + v)U:{p)] - $:T,expM^,(p;e,s)]. (2.8) 

n=l 

Here 

^pi/(P;0,-S) = Ss.ni^p{p) + Sin'&pi.P) + S2n^p{p) 

+S2ni>v{P + P) + S^n^uiP + P) + S^n^uiP + P) (2.9) 
-Sin^p{p + S5n^p{p + u) - S6n^p(p + l^) 

-S6n^u{p) - Syn'&uip) " Ssn'4'Ap)- 

The computation of the non-vanishing 2^ traces T„ shows that they all have the values 
±1/16. This structure is easily understood. Applying the projectors -Ps(3) and Ps(2) to 
a given vector successively leads to alternating projections on different subspaces. Each 
projection gives a factor ±l/-\/2, with the sign depending upon the There are 144 
values of n with T„ = 1/16 and 112 with T„ = —1/16. Obviously J2n Tn = 2 as it should. 
A further simplification is obtained from the symmetry T„(sj„) = T„(— Sj„), which is 
easily seen by inserting cticti between all projectors and using criPs{'m)ai = P_5(m) for 
m — 2,3. Together with the symmetry A'^^{p; 9, Sin) = —A'^^{p; 6, — Sj^) and the reality 
of the trace in (2.8) we can therefore fix, e.g. Sgn — 1, multiply (2.8) by 2 on the rhs, 
and restrict the sum over n to the 2^ values • • • , sy^. It is convenient to introduce 
sign factors e„ = ±1 and write 

r„ = ^ = rr[(P,,„(2)P,,„(3)) (P,3„(2)P,,„(3)) (P,,„(2)P,,„(3)) (P,,„(2)P,3„=i(3))]. 

(2.10) 

According to the previous remarks we dropped redundant projectors and fixed Ssn — 1- 
We may now write (2.8) as sum over the restricted set n in the form 

1 2^ 

TrUp^ip) = -Y^enCos[A;,{p;e,s)]. (2.11) 

° n=l 
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In order to perform the integrations over all the angles ip, "i?, (p we proceed essentially as 
in The various exponentials are expanded into a series of modified Bessel functions 
according to the formula 



exp[ezcosA]= ^ e^Ii{z)e^^^ for e = ±1. 



(2.12) 



l=—oo 



The replacement e ^ — e is obviously equivalent to the shift A ^ A + tt. Using (2.11) 
and (2.12) we thus find 



exp 



16, 



e„cos[v4;^(p;e,s)] 



n exp 



Ppu,n 



^e„cos[A;j,(p;e,s)] 



(2.13) 



= n E (e„)'/KY^)exp[.M^,(p;e,.)] 

= E n(en)'"^^'^%.(P)(|:)exp[^/;.(p)A;,(p;0,.)]. 

In the last step we exchanged the order of the product and the sum. The summation pa- 
rameters I must then be distinguished by indices referring to the corresponding plaquette 
Ppu and the index n; the l^^ip) run independently from — oo to oo. 

We can now write down the expectation value for a Wilson loop. We characterize it 
by a closed current loop J\{q) which is ±1 if the current runs in or against the direction of 
the link qx, and otherwise. For simplicity we exclude loops which run multiply through 
some links. For the calculation of the trace Tr[Vr(J)] we use again the decomposition 
(2.6) for the link variables on the loop, parametrized by (2.3). For every link qx on 
the loop we have two parameters s\{q) = ±1 (not three, because again neighboring 
projectors have to coincide), the whole set of these will be denoted by s. For a loop of 
lengt L we have 2L parameters s\{q) = ±1 and in total a sum of 2^^ terms. We count 
these by an index n, and denote the parameters by s\f^{q). The traces are called Wn. 
The Wn have a similar structure as the T„; their values are ±1/2^, and J2h = 2. The 
partition function now reads 



Z[J] = j EWf.eM^j:Uq)siM<dliq)] 

h qXb 

E ( n i^ny"^'^'^ii^Md) ^wK.ip)A;Ap; 0, s)]] 
n(mw)n^)- (2.14) 

Obviously the first sum over n is the expansion of the loop. The sum over lpi^{p) in 
the second line contains the action transformed as in (2.13), while the product over r,n 
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in the third hne contains the integrations together with the Haar measure. We may 
factorize the product over Ppi,,n and apply the addition theorem for the exponentials. 
This results in 



Z[J] = J EW^^M^EMq)sUQ)QxiQ)] 

h qXb 

nf^(w)n^^)- (2-15) 

We are now ready to perform the angular integrations over 0^(r), remembering the defi- 
nition of Apj^{p; 6, s) in (2.9). The integrations over ipfii^) and '^^{r) lead to a Kronecker- 
S which gives a constraint, while the 'i9^(r) integration involves the Haar measure (2.4) 
and leads to a more complicated function. We shall also call it a constraint for simplicity. 
It is convenient to define a symbol 6"'{C) for integer C by 



S"iC) = { p H-^QV^f^'?^ - / V(l - CV4) if C is a muhiple of 4, 1 . ^ _ . 

J-n ^ 2.-]q otherwise J ^ " ^ 

(2.16) 



The argument of the constraint which arises from the 0"(r)-integration becomes 



C;{r) ^ hsA,2)j;Ar) - Sie,,A)j;uir - u)] (2.17) 

y>fj.,n 

- E [s{8,7,6)nC^(r) - S(2,3,4)nC^(r " u)] + J^(r) 5;,^ (r ) , 

i/<fi,n 

where one has to use the first, second, or third subscript on s for a = ip,^,'^ respectively. 
So we end up with the following expression for the expectation value of the loop: 

z[j] = E ( n i^ny"'^'^iiu,)d)] nnc^ir)]. (2.18) 



3 Integrating out unnecessary variables 

The constraint equations (2.17), as they stand, have a quite different character than the 
corresponding ones in the abelian case in |jl|]. There one had one parameter l^v{r) for 
every plaquette r^^, in our case we have 2^ parameters l'^y{r) which are characterized by 
the additional index n. Any attempt of a physical interpretation of the l'^y{r) would be 
premature at this stage. 

Let us consider a fixed plaquette r^^ for the moment, and suppress the indices r^^. We 
first look for a suitable linear transformation from the parameters Z" to new parameters 
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m*. There are 8 combinations of the /" which play a special role, namely the 8 sums 
Sn=i "5m^" which appear in the constraints (2.17) (remember that sgn was fixed to 1). 
We will choose these 8 combinations as new variables m*, i = 1, • • • , 8, ehminate the first 
eight and keep the rest of the parameters as they are. The 8x8 matrix made up of 
the Sij with i, j = 1. ■ ■ ■ ,8 will be denoted by S. 

One may consider the as a set of 2^ vectors S„ = (si„, • • • , s^n, Ssn = 1) with 
components ±1. Up to now we did not specify the order in which these vectors S„ 
should be enumerated. It is now convenient to choose the S„ associated with the first 
eight values of n in such a way that en — 1 for n = 1, ■ ■ ■ , 8, with e„ the signs of the 
trace (2.10). This can be done in many ways. To be definite we give our choice in the 
appendix. Our criteria were a small but non- vanishing determinant of S (it is 128 for 
the S in A.l), and a structure as transparent as possible in some of the equations below. 
The transformation finally becomes 



128 8 128 



= ^SiJ"" ^Y^SijP + Y,SiJ'' for i = l,---,8, 

n=l j=l a=9 

m° = r for a = 9, •••,128. (3.1) 



In matrix form the transformation reads 



We have split m into an eight dimensional "constrained" vector rric, and a 120-dimensional 
"free" vector rrif. Here S is the 8x8 matrix defined before, with Sij — sij. T is an 
8 X 120 matrix, if we enumerate the columns from 9 to 128 we have TjQ, = Sj^- 

There is an important restriction which has to be imposed on the transformation. It 
is this restriction which will later on lead to the area law. If the run over all integers, 
the same will be true for the rrf with a = 9, • • • , 128, as seen from (3.1). But it is not 
true for the m* with i = 1, • • • , 8. From the inversion of the transformation in (3.2) one 
finds that only those appear, which fulfill the condition 

8 8 128 

F = Y,{S-%im^ - Y.iS-%i E 5-^" = integer for j = 1, • • • , 8. (3.3) 

i=l i=l a=9 

A computation of the matrix elements {S~^)ji shows that they are integer or half integer, 
and that the sums over the elements of any row are integer. Because the Sia are ±1, this 
implies that the second sum in (3.3) is automatically integer. Therefore the condition 
simplifies to J2i=i{S~^)jim^ = integer for j = 1, • • • , 8. This restriction finally becomes 
equivalent to 

m* — even for all i, or ~ odd for all i = 1, • • • , 8. (3.4) 
Introducing the transformation (3.2) into the product in (2.18) one obtains 
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En(^n)'"/p(^) = EE ( n^^^(-Srn«)(:§) I m^ar'^Im^Y,), (3.5) 

with P(m*,m") = Z^j('S'^^)jim* — Z]a(>S'"^T)jar«". The sum over m* only runs over the 
subset which fulfills (3.4). The for a = 9, • • • , 128 do not show up in the constraints, 
therefore the summations can be performed. To do this we introduce the integral repre- 
sentation for all the modified Bessel functions, thereby partially reversing the previous 
step of expanding the exponents. 

^'-^'(t^) = / TT exp[e— COS7 + i/7] for e = ±1. (3.6) 

The nf^ appear in the exponent now, and the summations can be performed with the 
help of the Poisson sum formula 



^exp 



= 27r^fc-E7i(^-'n-a 

ka .7 = 1 



27rA;. 



(3.7) 



All 7q are integrated over the interval — tt < 7q < vt, thus exactly one contributes in 
the sum on the rhs. Furthermore the 7,^ appear only as arguments of periodic cosines 
in the remaining part of the integrand. Therefore we may simply use the 5-functions to 



eliminate the 7a by putting 7^ 
with 



E5=i lj{S-^T)j^ for a = 9, • • • , 128. We thus are left 



En(^n)'"/p(^) = E/^---^e-p[^^(7) + ^ E iAs-'),m% (3.8) 

where we have introduced the function 

8 128 8 

^(7) = Ecos7,- + E e^cos[£lj{S-'na]. (3.9) 

j=i 0=9 j=i 

In order to fulfill the conditions (3.4) we put 

= 2m* + /, with nf integer, and / = or 1. (3.10) 

For our choice of the transformation, the /-dependence of (3.8) becomes particularly 
simple and only involves 73. 

We now reintroduce the suppressed indices Ppi, into (3.8), (3.10), and insert the result 
into the partition function (2.18). We find 



z[J]- E^^ E E n 



(hi drfs 
27r " ' 27r 



exp 



:§A(7) + 2zY: 7AS-%m^p) + ^7sfAp) 



(3.11) 



r/xa 
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The constraints (2.17) now depend upon m* and /. For the formulation it is convenient 
to extend f^y{r) to an antisymmetric matrix. Thus f^i,{r) = 0, 1 for /i < u, and f^y{r) = 
0,-1 for n > u. The constraints then become 



C,»(r) = 2 5: [m(«/'2) (r) - rf^'^ (r - z.)] - 2 ^ [m(«/-«) (r) - ^(^3'^) (r - u)] 

+ ^A,/^,(r) + J^(r)S;jr), (3.12) 

where denotes the left lattice derivative. The tensor f^uij) will be recognized as the 
Z2 field strength tensor. As in (2.17) one has to use the first, second, or third upper 
index of fh^y for a = ip,{}, ip. 

There is a symmetry relation in (3.11) which underlines the importance of the Z2 
tensor fpu{p)- Replace /3 — >• — /5 and substitute 7j — > 7j + tt for all j (the integrand is 
periodic). Using the definition of ^(7) in (3.9) and the fact that J2^=i{S~^T)ja = 1 for 
all a, one finds that ^(7) reverses sign, i.e. f3A{'y) stays invariant. The second term 
in the exponent of (3.11) changes by a multiple of 27ii, because Y^^=i{S~^)ji = 6si, and 
rhpyip) is integer. Finally the third term changes by mfpvip)- Iii this way one finds that 
the bracket (■ ■ ■) in (3.11) is even in (3 for fpuip) = 0, and odd in /? for fpuip) = 1- 

Up to this point all formulae were exact. It appears tempting now to proceed as 
follows in the continuum limit of large f3. If the function ^(7) has an isolated maximum, 
the integrals over 7^- are dominated by the region where ^(7) becomes maximal, and the 
integrations over the 7j can be extended to the interval from —00 to 00. A quadratic 
expansion around the maximum would then lead to gaussian integrals. If all the 
were equal to 1, we would indeed have a simple maximum at 7j = for all 7j. This 
would correspond to the situation in the abelian theory and to the replacement Ii{z) —>■ 
^z^-i ^y/2iTz which was used in [U (these authors used by mistake e~' instead of 
the correct e~' which has, however, no consequences there). 

In our case the different signs of the ea change the situation drastically. One finds 
that the function ^(7) assumes its maximal value of 16 if 71, 72, 73 are arbitrary, and 
7j = for j = 4, ■ ■ ■ , 8. Even a quadratic approximation in 7^ for j = 4, ■ ■ ■ , 8, and 
fixed 71,72,73 is not possible because the matrix of the second derivatives has two zero 
eigenvalues. For 71=72=73=0 even three eigenvalues vanish. So it would be necessary to 
go to a higher order in the expansion; but then the integrations are no longer gaussian 
and cannot be performend. This is the way in which non abelian gauge theory protects 
itself from being solved analytically! 

Nevertheless the present formalism will clearly show, how the area law for large loops 
arises. We keep (3 arbitrary, not necessarily large, and first solve the constraints. 



4 Solution of the constraints 

It is convenient to rewrite the -(^-constraints in the form 

S'[Ctir)]= E '^'[4A:.(r)]5c.w_4fc,w,o. (4.1) 
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This introduces additional sums over the ki^{r), and factors 5''[4A;^(r)]. The i^-constraints 
now also appear in form of a Kronecker-5. 

Let us first consider the constraints modulo 2, which obviously only concerns the 
second hne of (3.12). The factors sj^^(r) = ±1 may be dropped, and for all three cases 
a — ipj-djip we find the equations 



These are identical to the equations for /^jy(r) in the abelian case, except that they 
are equations modulo 2. They show already the appearance of a Z2 structure. We 
concentrate on a geometrical formulation of the solution. Recall that fnv{T) = 0, 1 for 
II < u and that ffj,u{r) is antisymmetric. Therefore it is convenient to use the symbol 
^ni^ — (1, —1, 0) for {/I < 1^, /I > ly, II — ly). Let jS" be a two dimensional surface, and 



The following statements hold: 

• If 5 has the Wilson loop as boundary, then fl£\r) is a solution of (4.2). 

• If is a closed surface, then fj£\r) is a solution of the homogeneous equation 
J2u ^ufj^^Hr) = (mod 2). The most general solution can be obtained as a 
superposition of a special solution fj^^ (r) with S bounded by the loop, and a sum 
over solutions of the homogeneous equation. 

The proof is obvious. Equation (4.2) involves exactly all the plaquettes which contain the 
link r^u. Links on the loop appear in an odd number of plaquettes of the associated 
surface S, while links which are not part of the loop appear in an even number 
(including 0) of plaquettes of S. 

In the following we will restrict the discussion to planar loops for simplicity. A 
solution of special importance is the layer belonging to the minimal surface of the Wilson 
loop W, 



2, if the loop is oriented appropriately. For different surfaces, on the other hand, this is 
not true. 

The general solution of the homogeneous equation J2u ^ufji^Hi") = (mod 2) can 
be written down explicitly. It depends upon the dimension d. In order to guarantee the 
antisymmetry of fnvir) we introduce the symbol (mod^^/ 2); it is identical with (mod 2) 
for // < I/, but reverses sign for // > i/. One then has 



5]A,/^,(r) + J^(r) = (mod 2) for all r,//. 



(4.2) 




e^t, if the plaquette r^^, is part of the surface 5", 
otherwise. 



(4.3) 





Ea ^nux^xfir) (mod^^ 2) for d = 3, 
xUir) (mod^^ 2) for d = 4. 



(4.5) 
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The function /(r) in three dimensions is unique up to a constant. For d = 4 one 
has a gauge freedom, i.e. adding a gradient AkA^t) to fnir) will not change fj^^^r). 
The simplest way to remove this ambiguity is to choose an axial gauge by imposing 
Z]k^k/k(^) = (mod 2). In both cases the values of /(r) and /^(r), respectively, are 

restricted to and 1. 

Switching from one surface S to another S' for the special solution can also be 
rephrased in terms of the solution of the homogeneous equation. In d = 3 dimensions 
it corresponds to changing /(r) by 1 inside the volume between the two surfaces. (The 
use of the left derivative in (4.2), (4.5) specifies which points of the surface have to be 
considered as inside or outside). For d = 4 one has to choose a three-dimensional volume 
spanned by the surfaces with, roughly speaking, normal vector in fi;-direction at the point 
r. One then has to change /^(r) by 1 inside the volume, and subsequently transform to 
the axial gauge. 

The essential part of the constraints has now been solved. We put 

Mr) = /(r"^(r) + Y: e,.x.AMr) {mod,, 2). (4.6) 

Xk 

The minimal surface layer fj^^^^'^r), defined in (4.4), is no longer a variable, but uniquely 
fixed by the loop. The /«(r) = 0, 1 are unconstrained. The index k on e^i^A/t, ^-nd 
in the sum appears for d = 4 only. For d — 3 it has to be dropped, here and wherever it 
appears in subsequent formulae. 

We next introduce the solution (4.6) into (3.12). The second line is now definitely 
even, therefore we denote it by 

'^RUir) ^ E^M^) + UrKnir). (4.7) 

V 

For later use we specify the variables upon which -R^ji(r) can depend. The term ^ju(r)s^^(r) 
is strictly local, i.e. only depends on the argument r. The /^(r), on the other hand, 
appear as A^AxfJj) with v ^ \. Therefore they enter also with shifted arguments r' . 
The points r' and r are neighbors in the sense that all components of r — r' are either 
or 1. Finally, one has to note that (4.6) is only an equation modulo 2. Therefore 

Ay] ,y{r) as well as -R^^i(r) can also depend on the minimal layer /^™"-'(r). 

The constraints (3.12) now become 



V>IJ, V<fJL 

+2Rl^{T)^45'^X{r). (4.8) 

For a = ip (first upper index i on m*) and a = ip (third upper index i on m*), i.e. for 
the even indices i, the index i appears in both sums of (4.8). For a fixed r one has 2d 
linear equations (corresponding to a — ip,(f and /i — 1, ■ ■ ■ ,d) ior 4d{d — l)/2 quantities 
m^^/'^'^^ . These equations are not independent due to the identity 

UCtir) - Clir - ^,)] = 2E[<(r) - - i^)]- (4.9) 

IJ, n 
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The rhs of (4.9) vanishes because /^jy is antisymmetric, and because the neighboring 
projectors in the loop have to coincide as mentioned in sect. 2. We checked exphcitly for 
d = 3,4 that the equations may be simply used to eliminate some of the Trf^^,{r). Any 
eliminated ?Ti^j,(r) depends linearly on other unconstrained m^'^i^^^r') and on RPf^{r'), 
where the components of r'-r are either or ±1. 

For a = d (second upper index i on m*), i.e. for the odd indices z, the situation is 
even simpler. Each index i enters only in one of the sums in (4.8), it is convenient to 
ehminate some of the rhj^^{r) and ml^{r). 

Finally we can write the solutions of the constraints in the following form, which 
eliminates some of the m^j,(r), leaving the rest unconstrained. 

r^,uir) = L;,[m;:,,,(r'),i?^,a(r')] +2A;;,(r), (4.10) 

with k^^yir) = (A;^(r), /cjy(r), 0) for {i = l,i = 7,otherwise). The L^^ are linear com- 
binations of their arguments, only coefficients 0, ±1 appear. The arguments r, r' are 
neighbors in the sense explained before. The loop current J^(r) enters only in i?^,^(r'). 

Note the drastic difference in the type of the constraints (4.2) (or the corresponding 
constraints in eq. (6) of ref for the U{1) case) on one hand, and the constraints (4.8) 
just considered on the other. The former involve a difference operator applied to one 
plaquette variable ffj,u{r), the corresponding Green function being non-local and couphng 
the solution to the current over a long range. In contrast, the latter constraints involve 
several plaquette variables m^f^^lr) and can just be used to eliminate some of these. This 
elimination leads to an almost local coupling to the current, involving neighbors only. 

5 Confinement 

The essential feature, which finally arose in our formulation, is the presence of the Z2 
field strength tensor ffj,u{r) which obeys the field equation (4.2). The solutions of this 
equation can be characterized by two-dimensional surfaces; a layer with the Wilson loop 
as boundary, possibly together with closed surfaces. One may expect that the presence of 
such a layer will lead to an area law. For a qualitative understanding of the confinement 
mechanism we use the explicit form (4.6) for the solution of the field equation (4.2). 
It contains the fixed layer /^™"''(r), together with the unconstrained Z2 variable fnif)- 
The solutions of the remaining constraint equations (4.8) for the fh^^^{r) have the form 
(4.10). 

Consider now the expression (3.11) for the partition function Z[J], and introduce the 
solutions (4.6), (4.10) of the constraints into the exponential on the rhs. This gives 

2E7,(5-'),.m;,(p)+78/p.(p) 

= 2 iAs-')Ai^UKAp')^ Rp'nip')] + ^K^ip)} (5.1) 

+ 78{/ir"nP) + Eep<^A.A,/.(p) (mod^, 2)}. 
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The Wilson loop enters into this expression in two different ways. First there is a 
dependence on the current Jp'{p')Sp,f^{p') which arises from the second term of Rp,f^{p') 
in (4.7). Secondly there is a dependence on the minimal layer which enters into 

the first term of as well as exphcitely in the factor of 73. 

The current J is present on a one- dimensional set, the minimal layer on a 

two-dimensional set. Besides this, both quantities enter in a quite similar way into (5.1). 
One may therefore expect, that for large loops the dependence on the one- dimensional 
current J can be neglected compared to the dependence on the two-dimensional layer 

y (mm) 

If we neglect the dependence on Jpi {p') the partition function becomes independent of 
the s'^pifiip') and the sum = 2 can be performed. Assuming that the 7- integrations 

have been done, the degrees of freedom are now in the remaining unconstrained rn^py{p) = 
—00, ■ ■ ■ ,00, the fnip) = 0, 1, and the A;^(r) = —00, ■ ■ ■ ,00 introduced at the beginning 
of sect. 4. 

The discussion of (4.8) showed that the solutions couple neighbors only. This means 
that (5.1), which appears in the exponential in (3.11), only depends on these variables 
with arguments p,p',p"', here p',p" are neighbors in the sense that all components of p'-p 
are 0, ±1, all components of p"-p are 0, ±1, ±2. 

We digress for a technical point. Neither the exponential in (3.11) with it's com- 
plex argument, nor the factors 5''[4fc] are positive definite. Actually, according to the 
definition (2.16), one has Z]fc'^''[4/c] = 0, because the Haar measure fulfills i^(0) = 0. 
If desired, one could bring the expression into the usual form of a partition function 
with positive summands, by performing a twofold partial summation with respect to the 

The whole loop dependence is now in the belonging to the minimal surface. 

It acts like a space-time dependent external field, comparable, say, to a constant magnetic 
field switched on in a finite volume of an Ising model. Z[J] is a partition function where 
the variables couple to neighbors only, a well known standard situation in statistical 
mechanics. For large subsystems it therefore factorizes into products refering to the 
subsystems and, correspondingly, has an exponential dependence on the volume. This 
is, of course, nothing else but the fact that the free energy is an extensive quantity. 
Rigorous proofs, which apply for any dimension, can be found in 0. 

Consider now a loop < a;i < i?, < X4 < T in the Xi-X4-plane for definiteness, 
with R and T large. We divide the Xi-X4-plane inside, as well as outside of the loop, into 
rectangles; these rectangles are then extended to c?-dimensional boxes into the orthogonal 
directions. This means that we define regions V^'-"^ by the inequalities rj"'' < xi < 
r2"\ti"^ < X4 < t'2'\x2,x^ arbitrary. The rectangle r^"-* < Xi < r2^\ti^^ < X4 < 4"'' 
has to lie either completely inside the loop, or completely outside the loop. If not only 
R and T, but also all the differences r2^^ — r^"'' and ^2"'' ~ 4"^ are large, the partition 
function will factorize, 

Z[J] = llZ^''\ (5.2) 

n 

Consider now the ratio Z[J]/Z[0], with Z[0] the expression without loop. Obviously all 
the outer factors cancel. For the inner ones, on the other hand, one has /p™"Hp) — 1 
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in the numerators, but in the denominators. Thus the ratios are different from 1. 
Because of the factorization property, the volumes of the regions F*^") have to enter in 
the exponent. This finally implies that the area A — R x T of the loop enters in the 
exponent, so the result may be written as 

Z[J]/Z[0]=exp[-aA]. (5.3) 

Wc have obtained the law for large loops. 
Several comments are appropriate here. 

First one may wonder what would happen with our argumentation, if we would 
replace /^^'"^(p), associated with the minimal surface, by a solution f^^\p), belonging 
to a different surface S. Obviously the simplicity of the situation for the regions 
inside and outside would break down, and factorization would not lead to a simple 
relation. The minimal surface is really unique for the argumentation. 

The neglection of the dependence on J would certainly have been wrong if performed 
in the original expression described by the Euler angles %l)^{r),'&^{r), ^n{r). There one 
had only the current J but no layer showed up. Therefore confinement has to 

arise from the dependence on J in a complicated way. In our formulation the formalism 
led to another quantity, the minimal layer This appears as the natiu^al quantity 

which describes the long distance physics and dominates the residual direct dependence 
on J. 

For illustration one can have a look on the strong coupling limit. According to the 
discussion at the end of sect. 3, the bracket (■ • ■) in (3.11) is even in j3 for fpu{p) = 0, and 
odd in f3 for fpu{p) = ^pu- Therefore the order only contributes outside the surface 
layer, while the order (3 terms come from plaquettes on the surface layer. In this way 
we recover the well known lowest order strong coupling result Z[J]/Z[Qi\ ~ More 
important, we have seen that indeed is the crucial quantity. 

For a Wilson loop in the adjoint representation one does not expect an area law, 
because the charges can be screened by pair creation. This can be easily checked in our 
approach. The traces in the adjoint and in the fundamental representation are related 
by TrW(^i) = [TrVF(i/2)]^ — 1. With our parametrization we obtain 



T 2 



qXb 



WnWn> exp 



qXb 



(5.4) 



The last term on the rhs of the modified equation (3.12) becomes Jtj,{r)[s1j^f^{q) + ^^^^/(g)] 
and is always even. Therefore (4.2) becomes a homogeneous equation, no minimal layer 
and no area law will appear. Similarly one can see that we don't get confinement if we 
replace SU{2) by SO (3). 

With some technical effort or a more streamlined approach it should be possible to 
carry through a similar analysis for SU{3). It would be nice to see, how the formalism 
would create the expected Z^ structure. 
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Our conclusions which led to the area law would break down if the result, by some 
reason whatsoever, would be independent of fpu^p), thereby giving a vanishing string 
tension. This appears hardly possible for a finite lattice. We have seen before that 
there is indeed an essential dependence on fpuip) in the strong coupling limit /? — 0. 
Such a dependence must survive for all finite (3 because the original expression Z[J] in 
(2.1) clearly fulfills the strict inequalities < Z[J] < Z[0]. The string tension might, 
however, vanish in a certain region of (3 after performing the thermodynamic limit. In 
particular such an effect could be expected in higher dimensions, where the presence of 
the two dimensional layer becomes relatively less important than in lower dimensions. 
Indeed it is known that lattice SU (2) has a first order phase transition for d = 5 at 
13c = 1.642 ±0.015. 

We come back to d = 4. At the end one is interested in the continuum limit /? — > oo 
which requires a particular investigation. If the string tension is a physical quantity and 
f3 goes to infinity as prescribed by the renormalization group, a non vanishing string 
tension for the lattice theory will persist in the continuum limit. 

6 Interpretation and conclusions 

There is an extensive literature on the various pictures of confinement which cannot be 
discussed here. For a recent review we refer to 0. We come directly to the physical 
interpretation of our results. The key is equation (4.2) for the Z2 field strength tensor, 

Y,AJ^,{r) + J^{r)=0 (mod 2). (6.1) 

V 

The solutions in form of layers on two-dimensional surfaces were discussed in detail in 
sect. 4. 

In (i=3 dimensions put fpu{r) = J2x^fiu\B\{r) (mod^j/ 2). Then (6.1) becomes 
V X B = J (mod 2). The magnetic field B has sources corresponding to magnetic 
monopoles. It is reasonable to use the right derivative in the divergence, and to associate 
the monopole density p with cubes as usual. We therefore define 

p(ri23) = E^l"'''*^^-W (mod 2). (6.2) 

The solution fj^^'^^^r) then immediately leads to a double layer of monopoles in the 
cubes on both sides of the minimal surface. 

For rf=4 we define the dual tensor /^^.(r) = {l/2)J2x,^(^puXKfxK{r) (mod^^, 2). The 
conserved (mod 2) magnetic current lives on 3-dimensional cubes r^xn, where p. A, k 
denote the three directions orthogonal to fi. 

Mrpx.) = E ^'u'''"'^ UAr) (mod 2). (6.3) 

V 

Consider a loop in the Xi-X4-plane, with 0:4 interpreted as euclidean time. Let 0:1,0:4 be 
within the loop and suppress the X4-extension of the cubes. For the solution Z^™"-* (r) we 
then get a non-vanishing on all plaquettes in the Xi-X2-plane and in the Xi-Xs-plane 
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which contact the hne X2=X3=0. We thus have a string of electric field Ei{r) = fuir) in 
Xi-direction, concentrated on X2=X3=0. This is surrounded by magnetic current loops 
parallel to the X2-X3-plane which circle around the electric string. The configuration 
is therefore just dual to an Abrikosov vortex in a type II superconductor, where the 
magnetic field is encircled by the electric current. Flux quantization is evident, the 
Z2-structure only allows for one unit of fiux. 

Configurations f^u{f) in (4-6) with /^(r) 7^ belong to other surfaces which are 
bounded by the loop as discussed in sect. 4. In addition closed surfaces can appear. 
The interpretation is similar as above. For illustration, connect e.g. two points in the 
a;i-a;2-plane by a path in form of a stair. Then lies on the plaquettes which point 
from the stair into positive and negative xs-direction. All surfaces are summed with 
the appropriate weight in the partition function. A careful investigation of the various 
weights should give information about the extension of the electric fiux tube. 

We finally check the dual London equation, V x J = -i-E (mod 2). It is easily seen 
that (V X J)^(r) (mod 2) is equal to (1) if the link has contact to an even (odd) 
number of plaquettes with no n- vanishing magnetic current. For the examples discussed 
above this means that V x J runs along the boundary of the set of plaquettes which 
carry the magnetic current. For the minimal layer, V x J is parallel to E as it should. 
It is, however, not concentrated on X2=a;3=0 as the electric field, but on the four lines 
X2=0, X3 = ±1 and X3=0, X2 = ±1- For the stair, V x J is shifted by X3 = ±1 with 
respect to E. In general the dual London equation is essentially fulfilled, the two sides 
of the equation are just slightly shifted against each other. This might be interpreted by 
a non-vanishing Ginzburg-Landau coherence length ^ which leads to a "normal" region 
near the string, where the London equation is not valid. 

Let us compare with some familiar types of monopoles in the literature. 

The charges of the U{1) monopoles in |l|] can take all integers, in obvious contrast to 
our Z2 structure which only allows and 1. 

A popular definition of monopoles in SU{2) is discussed e.g. in ^j. Let ?7(p) = 
sign Tr U(^p) denote the sign of the plaquette action, and rjc = Yipedc V{p) the product of 
the ?7(p) around the boundary of the cube c. Then r^c = — 1 represents a monopole in the 
(space like) cube c. There is a Z2 structure as in our case. 

Another frequently applied definition, reviewed e.g. in 0, uses the maximal abelian 
gauge. In a first step one maximizes the quantity R = J2rtiTr[a3Uf^{r)a3U^{r)]. The 
link matrices are then decomposed into a non-abelian and a f/(l) part, e.g. one can take 
the abelian link angle as the phase of [Ufj_{r)]ii. The U{1) monopoles are then defined 
according to the DeGrand Toussaint construction |TD| which allows monopole charges 
0,±1,±2. 

The monopoles which naturally arose in the present work have no direct relation to 
any of these. An unconventional feature of our approach is the presence of monopoles in 
every configuration. In the approaches mentioned above there are plenty of configurations 
without any monopoles, namely those near the perturbative vacuum. In our case there 
is always a surface S bounded by the loop. This is associated with an electric string and 
accompanied by monopole vortices. From a physical point of view this appears quite 
attractive. 



16 



Acknowledgement: I thank I. Bender and H. J. Rothe for reading the manuscript and 
for valuable criticism. 



A Appendix 



For definiteness we give here the matrix S used by us in sect. 3 when selecting a 
convenient subset of the S„. It reads 



S 



-1 
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-1 


1 


-1 


1 


1 


1 


-1 


-1 


1 


-1 


1 


-1 


1 


1 


-1 
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1 
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1 
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-1 


-1 


1 


1 
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-1 


1 


1 


-1 
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-1 


-1 


-1 
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-1 
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-1 


1 


-1 


-1 


1 


-1 


-1 


1 


-1 


1 


1 


1 


1 


1 


1 


1 


1 


1 



(A.i: 



Recall that the columns of S consist of 8 of the vectors S^, with the property that 
en = +1. The last component, corresponding to Sg, was fixed to 1. 
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